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QCD. In accordance with these ideas, we assume that the soft Pomeron intercept is rather large, and the 
<0 slope of the Pomeron trajectory is equal to zero. We derive analytical formulae that sum both enhanced 
T 7^ and semi-enhanced diagrams for elastic and diffractive amplitudes. We fit the available experimental data, 

and predict the value for cross sections at the energies accessible at the LHC. The main corrections to the 
^ model are studied and evaluated. 
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1. Introduction 

Thanks to the hard work of experimentalists and theoreticians over the past four decades, we know that the 
most economical and reliable method for describing soft interactions at high energy, is the phenomenology 
based on the soft Pomeron and secondary Reggeons (see Refs. [1-3] for details ). Consequently, we believe 
that the future theory should be a theory of Pomeron and Reggeons and their interactions. However, 
numerous attempts [2, 4] to build such a theory have failed, as one was not able to specify the Pomeron 
interaction, as well as interaction of Pomerons with the target. We are doomed to make ad hoc assumptions 
about the vertices of multi- Pomeron interactions (see Ref. [4-6,9,10]) which specify the approach, but do 
not make it more theoretically reliable. Such assumptions seem unavoidable, as there is no theoretical 
approach to non-perturbative QCD and, on the other hand, soft high energy processes appear to be 
typical examples of non-perturbative physics at long distances. These processes are even more difficult 
to calculate, as approximate methods such as QCD sum rules and/or effective theories, as well as the 
lattice QCD approach, cannot be employed to determine the high energy amplitudes. The success of 
any phenomenology does not establish a theory, however over the past two years a new approach has 
been developed (N=4 SYM), which allows one to study theoretically the regime of the strong coupling 
constant [11]. 

1.1 N=4 SYM 

We use this theory as a guide for the physics phenomena occurring in this regime. The attractive feature 
of this theory is that N=4 SYM with small coupling, leads to normal QCD like physics (see Refs. [12,13]), 
with OPE and linear equations for DIS, as well as the BFKL equation for the high energy amplitude. The 
high energy amplitude reaches the unitarity limit: black disc regime, in which half of the cross section 
stems from elastic scattering, and half relates to processes of multiparticle production. 
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On the other hand, making use of the fact 
that with the aid of the AdS/CFT correspon- 
dence [11] this theory can be solved analytically, 
and can be reduced to weak gravity in AdS§ 
space. 

In the strong coupling limit the following 
main features of this theory citeBST,HIM,COCO,BEPI,LMKS 
are manifest: (i) it has a soft Pomeron, which 
in this case is the reggeized graviton with the 
large intercept a(0)jp = 2 — 2/vA where A = 
AirN c c?g M and a™ 7 is the QCD-like coupling; 
(ii) the main contribution to the total cross sec- 
tion at high energy is due to the processes of elas- 
pdsiii scattering and diffractive dissociation; (iii) 



Figure 1: The behaviour of the Pomeron trajectory in N=4 the leading Pomeron trajectory has a form shown 

B om4 
, namely, a Regge pole with otjp = in 

the scattering region (t < 0) while a'jp > in 

the resonance region of positive t; and (iv) the Pomerons (gravitons) interact by means of the triple 

Pomeron vertex which is small (at least oc 2/vA) • The minute value of a'jp is not related to the small 

t Iftpm4 

size of the partons in this theory. It is related to small values of the fifth coordinate r (see Fig. |i|) . The 
physical meaning of this coordinate is, related to the typical size of the colliding particles. It should be 
stressed that all these features are an integral part of the theory and therefore, for the first time we have 
theoretical justification for using Reggeon-type phenomenology in high energy scattering. 

In this paper we present our approach based on two major assumptions: it reproduces the main features 
of N=4 SYM, and it provides a natural matching with the perturbative QCD approach. For the sake of 
completeness we discuss the perturbative approach below. 



1.2 Matching with perturbative QCD 

In perturbative QCD the high energy amplitude has been calculated in the leading log approximation, 
in which < 1 but as logs ~ 1, where -y/s = W denotes the energy in the cm. frame [22-24]. This 
amplitude can be written as the exchange of the BFKL Pomeron [23], which has the following form for the 
cross section for the scattering of one colourless dipole of size r on another dipole of size R 



a (r, R;s) = J d 2 bd 2 p\ d 2 p2 



+ioo+e j 

— V{r, Pl -v) e^ Y V{R, P2 ;-y) S^( Pl -p 2 -b) (1.1) 



ioo+e 



2m 



where 



12 



uj(v) = a 5 (2^(1) - ^(1/2 + iv) - ipil/2-w)) 



1.2 



13 
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and ip(z) = dlogT(z), and T(z) is the Euler gamma function. 

One can see that the BFKL Pomeron is not a pole in angular momentum, but it is a cut. However, 
at fixed v which is the conjugate variable to In (r 2 /i? 2 ), it is a pole. The position of this pole does not 
depend on momentum transfer (or b) and, therefore, the corresponding a' is equal to zero. 

It turns out that in the wide range of energy l/a 2 s > Ins > 1 the scattering amplitude can be 
expressed as the sum of BFKL Pomeron exchanges and their interactions (see Refs. [25-32]. Perturbative 
QCD specifies the Pomeron vertices, and in the LO approach, the only vertex that contributes is the triple 
BFKL Pomeron vertex [31,32]. 



1.3 Models of high energy soft interaction 





We are incapable of building a theory 
of high energy scattering, as the problem 
of the confinement of quarks and gluons in 
QCD has not been solved, and we do not 
have a theoretical tool for describing the in- 
teractions of quarks and gluons at long dis- 
tances. Therefore, we are doomed to build 
models that take into account our ideas of 
the b ehaviour of QCD at long distances. 
I^eelieve that such models should include 



Figure 2: The exact Green's function of the Pomeron as a sum of everything that we know about QCD at 

enhanced diagrams and the exact triple Pomeron vertex. . . . . . . . 111 

long distances and, m particular, should 

absorb all that we have learned about high energy scattering in N=4 SYM. We can summarize this knowl- 
edge as a list of criteria that a model should satisfy, namely, 



1. 



The model should be built using Pomerons and Reggeons as the main ingredients; 

The intercept of the Pomeron should be rather large. In N=4 SYM we expect Ap = ajp(0) — 1 = 
1 - 2/VA « 0.3^-0.4, since the estimate for A from the cross section for multiparticle production as 
well as from DIS at HERA [33] is A = 5-^9; 



3. 
4. 



(0) 



0; 



A large contribution should come from processes of diffraction dissociation, since in N=4 SYM at 
large A only these processes contribute to the scattering amplitude. In other words, the model should 
include the Good- Walker mechanism [45] as the main source of the diffraction production; 

5. The Pomeron self-interaction should be small (of the order of 2/\f\ in N=4 SYM), and much 
smaller than the vertex of interaction of the Pomeron with a hadron, which is of the order of A; 

6. The last requirement follows not from N=4 SYM, but from the natural matching with perturbative 
QCD: where the only vertex that contributes is the triple Pomeron vertex. 
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Figure 3: Examples of 'fan' Pomeron diagrams: in Reggeon Calculus (the first three diagrams) and in the generating 
function approach (the last three). 

In this paper we continue to develop a model that satisfies all above criteria. In our previous paper 
(see Ref [7]) we suggested a model that includes the large diffraction from Good- Walker mechanism and 
a Green's function of the Pomeron that sums all enhanced diagrams (see Fig. However, we neglected 
the contribution of other Pomeron diagrams (see for example the so called fan diagrams in Fig. Hfj^The 
line of argument for justifying this was the following: The value of the triple Pomeron vertex turns out 
to be small, and the cross section for the diffraction production in the region of large mass makes up 10 - 
20% of total cross section of the diffractive production. Therefore, we decided in the first approximation 
to neglect the process of diffraction of large masses, and calculate them perturbatively. In this paper we 
sum all 'fan' diagrams together with the enhanced ones. 



2. Main ideas and assumptions 
2.1 Soft Pomeron 

As we have mentioned we wish to suggest a procedure which contains the main features of N=4 SYM, 
and provides a transparent matching with perturbative QCD, which plays the role of the correspondence 
principle in our approach. 



The good news is that the soft Pomeron is a natural ingredient in N=4 SYM. Actually, in N=4 SYM 

mpm4 

es (see Fig. Ill) . Our simplification is to replace this set of 
m red. We feel confident with the assumption regarding 



we have an infinite series of different Regge po 
poles by one pole: the Pomeron, shown in Fig. 



the existence of the soft Pomeron, since it has been successfully utilized in high energy phenomenology for 
the last forty years, and it appears in QCD (see above). 

^From N=4 SYM we expect that the value of A = ajp(0) — 1 = 1 — 2/y/X could be large. 
2.2 ol'jp= 

Our second ingredient is that the slope of the Pomeron a'jp = 0. As we have discussed, a'jp = is in perfect 
agreement with N=4 SYM predictions. It also agrees with the recent fit to high energy data [7] . This 
fact might appear strange, as for a long time the widely accepted value for a'jp was a'jp = 0.25 GeV~ 2 [43]. 
Comparing with the parametrization of Ref. [43], two new ingredients have been introduced: a large value 



- 5 - 



of A ~ 0.3 instead of A ~ 0.08 in Ref. [43], and sufficiently large shadowing corrections, which were 
considered small in Ref. [43]. 

We believe that a small value of a'p, is a signal that rather short distances contribute to the soft 
interaction at high energy, in agreement with all approaches to the Pomeron structure considered above. 

It should be stressed that N=4 SYM gives a new possibility for high energy assymptotic behaviour: a 
Regge cut (Pomeron) that does not move in the scattering region, while it has a slope a'jp in the region 
of positive t. In Ref. [44] such a possibility was missed. It is interesting to note that in N=4 SYM this 
property of the Pomeron (graviton) contribution, stems from the integration over the fifth coordinate z, 
which has the meaning of integration over the possible size of the hadrons. 



2.3 Good- Walker mechanism 



The third conclusion that we derive from N=4 SYM, is the large contribution of the diffractive dissociation 
processes. In our approach the diffraction dissociation is taken into account using the two channel model 
which we have developed in a number of papers (see Ref. [7] and references therein). In this formalism, 
diffractively produced hadrons at a given vertex are considered as a single hadronic state described by the 
wave function typ, which is orthonormal to the wave function ^ of the incoming hadron (proton in the 
case of interest), < ^h\^D >= 0. We introduce two wave functions ipi and ip2 which diagonalize the 2x2 
interaction matrix T 

Ai t k =< 4>i ipkl^lA' tpk' >= A,k Si,i' h,k'- (2-1) 

In this representation the observed states are written in the form 



where, a 2 + j3 2 = 1. Using Eq. 



ifj h = a ipi + j3 V>2 , 

i>D = -0 ipi + a 1P2 , 
J, we can rewrite the unitarity constraints in the form 

ImA i>k {s,b) = \A ijk (s,6)| 2 + G™ k {s,b), 



M 



(2.2) 
(2.3) 

(2.4) 



where Gft is the contribution of all non diffractive inelastic processes, i.e. it is the summed probability for 

these final states to be produced in the scattering of particle i off particle k. 

IJHX T 

A simple solution of Eq. (|2.4[) has the same structure as in the single channel formalism, 



1 — exp 



(2.5) 



G™ k (s,b) = l-e W (-n hk (s,b)). (2.6) 

^From Eq. (p we deduce, that the probability that the initial projectiles (i, k) reach the final state 
interaction unchanged, regardless of the initial state rescatterings, is Pf k = exp (—Q i k (s,b)). 

The opacities Qi k (s,b) have to be determined and we will discuss this below. In general this two 
channel approach is a particular case of the Good- Walker mechanism [45] for diffractive production, and 



2CHM 



2CHM31 



2CHM32 



UNIT 



2CHM1 



2CHM2 
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we can account for diffraction in the region of small mass using this approach. For the region of large 
values of produced mass in diffractive dissociation, it is necessary to develop a theoretical approach to the 
Pomeron-Pomeron interaction. 



2.4 Pomeron -Pomeron vertices 

The general framework for accounting for Pomeron interactions was developed in 70's (see Refs. [2,46] in 
the framework of Reggeon Calculus. However, being purely phenomenological, Reggeon Calculus is not 
able to fix the vertices of Pomeron interactions. The values of vertices, and their number, have to be 
introduced in Reggeon Calculus from the microscopic theory. 

Our third ingredient is the assumption that only the triple Pomeron vertex is essential. This assumption 
is in full agreement with N=4 SYM and perturbative QCD: in both these theoretical approaches only this 
vertex contributes. We consider this to be essential, as this assumption provides a natural bridge to 
perturbative QCD. To illustrate this matching we consider the sum of the 'fan' Pomeron diagrams (see 

jfanst 

Fig.p^ 

Instead of the generating functional for high density QCD (see Ref. [36] ), we can introduce the 
generating function: 

oo 

Z(Y-Y ;u) = ^2P n (Y-Y )u n (2.7) 

n=l 

where P n (Y — Yq) is the probability to find n-Pomerons at rapidity Y — Yq. By fixing the size of the 
interacting dipoles, one can see that Eq. differs from the generating functional in the dipole approach 
of high density QCD. 

The full set of the 'fan' diagrams can be summed using the following equation for the generating 
function 

9Z %Y Y °' ,U) = - An ( 1 " u ) Z ( y " y o;^) (2.8) [zfeq; 

with the amplitude that is given by the following equation 



oo . n Aj oo 

N (Y; H}) = " E / ^( y o) n^^(^. «) l«=l = " £(-l) n 7nCro) P(Y - Y ) (2.£ 



n=l i=l n=l 



ZF 



NF 



If we neglect dependence on the size of the dipoles, then both equations coincide with the equation in the 
dipole approach to QCD. These equations sum the set of 'fan' diagrams, if we replace (see Fig. Hf 5 

91 -+ h = G 3]P ->■ A ; g 2 -> 7 = (2.10) [rep] 

where A denotes the intercept of the Pomeron, and 7 the amplitude of the interaction of the 'wee' parton 
(colorless dipole) with the target. The scattering amplitude is equal to g\ N (Eq. (pT9|)). 

Assuming that only the triple Pomeron vertex is essential, we achieve the matching between pertur- 
bative QCD and our approach. We will show below that this matching can be demonstrated for sets of 
Pomeron diagrams, more complicated than the 'fan' diagrams. 
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The theory that includes all the ingredients that have been discussed above, can be formulated in a 
functional integral form [31], 

Z[$,$+] = J D<$>D$ + e s with S = S + Si + S E , (2.11) 

where So describes the free Pomerons, Si corresponds to their mutual interaction, and Se relates to the 
interaction with the external sources (target and projectile). Since a'jp = 0, So has the form 

So = J dY$ + (Y) j- A + aU(Y). (2.12) 

Si includes only triple Pomeron interactions and is of the form 

Si = 93ip j dY ($(y)$+(y)$+(y) + h.c.) (2.13) 

Se depends on our model for the interaction of the Pomeron with the scattering particles, which we will 
specify later. 

This theory, as any theory of Pomeron interactions, is written in such a way that the high energy 
amplitudes satisfy t-channel unitarity. However, s-channel unitarity remains a problem. In Refs. [37] it 
was shown that to satisfy s-channel unitarity we need to add to the interaction term (Si), the four Pomeron 
vertex ~r(2 -> 1)($ + ) 2 ($) 2 . 

For a better understanding of s-channel unitarity we reformulate the theory, given by the functional 
integral of Eq. ( jzTLl ), in terms of the evolution equations for the system of partons. As we have mentioned, 



SO 



SI 



for perturbative QCD these partons are colourless dipoles, as was shown in Ref. [34], that can decay and 
merge: one parton to two partons and two partons into one parton, with probabilities T(l — > 2) and 
r(2 —> 1) respectively. For such a system of partons, we can write a simple evolution equation (Fokker- 
Planck equation). Indeed, let P n (y) be the probability to find n-parton (dipoles) with rapidity y in the 
wave function of the fastest (parent) parton (dipole), moving with rapidity Y > y. For P n (y), we write 
down a recurrence equation (see Refs. [36,48]) 

8 P (v) 

_ nyyi = r(1 ^ 2) { _ nPn + ( n _ 1)Pn _ l} + r(2-H) {-n{n-l)P n + (n + 1) n P n+1 } . 

oy 

(2.14) PNEQ 

In each bracket the first term on the r.h.s., can be viewed as a probability of a dipole annihilation in the 
rapidity range (y to y — dy) (death term). The second is a probability to create one extra dipole (birth 
term). Note the negative sign in front of dP n (y) / dy. It appears due to our choice of the rapidity evolution, 
which starts at the largest rapidity y = Y, of the fastest dipole and then decreases. The first two terms are 
responsible for the process of parton decay, while the last two terms describe the contribution of partons 
merging. 

I 7F I PMF.n 

Using Eq. ( |2. y[ ) we can re- write Eq. Q2.14 ) in the form 

dZ(y,u) , . dZ(y.u) , . d 2 Z(y, u) . , , , 

= -ku(1-u) + u(l -u) ^ L - L . 2.15 [GFEQ^ 

oy ou a l u 1 1 
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where y = Y{2 ->■ 1) {Y - y) = A 7 (Y - y) and k = I/7. 

Eq. (^^) should be added with the initial condition at Y = Yq, as well as the boundary condition 

Z(y, u = 1) = 1, (2.16) 



which follows from the physical meaning of P n as a probability. It turns out that Eq. (12.151) can be solved 
analytically with an arbitrary initial condition [49,50]. Therefore, formally speaking we could derive the 
initial condition for different processes and obtain the solution using the approach of Ref. [50]. However, 
we chose a different and more traditional way to tackle the problem. As was suggested in Ref. [46] the 
solution of the Pomeron interaction problem should have two steps. Step one, is to find the correct Green's 
function of the Pomeron, and the exact vertex of the triple Pomeron interaction, by summing all enhanced 
diagrams (see Fig. ^ Step two: to solve the problem of the interaction of new (exact) Pomerons. One 
expects that the interaction in terms of the exact Green's function and vertices, will be much simpler at 
high energy, as the main effect of the interaction has been taken into account in step one. 

It is necessary to require that there be only one fastest parton (dipole), which is Pi(y = Y) = 1, 
while P n> i(y = Y) = 0. In this case we have the following initial condition for the generating function 

Z(y = Y) = u. (2.17) 

and the solution to Eq. (pTTo) will give the exact Green's function of the Pomeron. However, we prefer to 



develop a different technique for finding the Pomeron Green's function, which makes the calculation more 
transparent, and leads to explicit analytical formulas for physical observables (see section 3.1). 

2.5 The phenomenological parameters and their typical values 

Unfortunately, even with all assumptions that we have made, our approach is still phenomenological, since 
we have to determine the parameters of our interaction from a fit to the experimental data. However, 
using our main idea that the soft interaction stems from rather short distances, we are able to give some 
estimates for these parameters. 

First, we list all of these parameters: 

1. For description of Good- Walker mechanism of diffraction production in the two channel model, we 
need two phenomenological functions ffi(6 ) m h! g2(b), whi ch describe the vertices of interaction^ of 



the Pomeron with state 1 and 2 (see Eq. (2T2)^md Eq. ( 12.20 ), and one number (3 see Eq. (|2,zj) an d 



Eq. (g^TTr 

2. The Pomeron intercept A = a(0) — 1 ; 

3. The low energy amplitude of dipole -target interaction 7; 

Since we believe that the short distances contribute to the Pomeron structure we expect that 

„2 . 



INC2 



INC1 



A oc as; 7 oc a s ; g\ k> g 2 oc 1 (2.18) | EST 
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For gi(b) we use the entire phenomenological assumption 

9i(b) = giS{b) 



(2.19) [S 



where 5(6) is the Fourier transform of the dipole formula for the form factor 1/(1 + q 2 /m 2 ) 2 . 



3. Summing the enhanced diagram 

In this section we sum the enhanced diagrams, and obtain the exact Green's function of the Pmeron, as 
well as the exact vertex for the triple Pomeron interaction. To achieve this we employ the approximation 
that has been developed in Ref. [53] (MPSI approximation), and the improved version as in Refs. [54,55], 
so as to adjust this method for the summation of the Pomeron loop diagrams, for Pomerons with the 
intercept A > 0. 



3.1 Improved MPSI approximation 

To illustrate the method, we calculate the first enhanced diagram of Fig. 

*Y ryi 



hi 



A (Fig. |T= -g l92 Gi F / dy x I dy 2 G(Y - Vl ) G z ( Vl - y 2 ) G(y 2 - 0) 

Jo Jo 

-2 / ^. I d A(Y +yi -y 2 ) 



= ~ 9i92G 31P I dyx I dy 2 
Jo Jo 

~ ~ [22AY , 2 Ay , a 2v AVl 

where G^jp = A 7 and gt = gij y/7 (see Fig. ^ffor the notation) 



A 2 



(3.1) 





Y-yi~ 1/A 





- V2 - « 1/A 

Figure 4: The first enhanced diagram for the Figure 5: The first enhanced diagram in the form 



enhl 



Pomeron with intercept A > 0. 



enhlgf 



suited for the MPSI approximation. 



II 



The main idea of the MPSI approximation is to take into account only the first term in Eq. ( |3.l| ) 
neglecting other terms, since they are suppressed as exp[— Ay]. This term is the result of integration for 
Y — yi ~ 1/A and y 2 — ~ 1/A (see Fig. Hp l The general expression for the sum of the enhanced diagrams 
in MPSI approximation is shown in Fig. 
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One can see that the MPSI approximation is the 
i-channal unitarity constraint adjusted to Reggeon Cal- 
culus, in the form of the generating functional ( gener- 
ating function in our case for summing Pomeron inter- 
actions). From this picture, and from our knowledge 
of the cascade described by the 'fan' diagrams, we can 
find the answer for the sum of enhanced diagrams only. 
The physical meaning of the introduced parameters is 
clear: 7 is the low energy amplitude for two partons 
(dipoles) scattering at an arbitrary rapidity Y', and A 
is the value of the vertex for the decay of one parton 
(dipole) to two parton (dipoles). It should be stressed 
that the answer does not depend on the value of Y', 
but it should be chosen somewhere in the central region 
gfSPHife scattering. 



Figure 6: The exact Green's function of the Pomeron 
as a sum of enhanced diagrams in perturbative QCD 
in the MPSI approximation. 



To find the generating function that describes the 

JZEEQ 

sum of the 'fan' diagrams we need to solve Eq. 



J- 



This is not a difficult task, and for illustrative purposes 
we will do it in a different way. First we introduce a new generating function 



N(Y-Y',y) = 1 - Z(Y-Y',u = 1 - 7) 



(3.2) 



MPSI2 



^From Eq. ( fO| ) one can associate ./V with the scattering amplitude, when the variable 7 is equal to the 
low energy parton amplitude. If we introduce a new variable 



1R 



7 



1-7 



(3.3) 



Eq. (|2.6[) reduces to the form 



dN{Y-Y'; lR ) dN(Y-Y'- lR ) 
= — 



d (Y - Y') 



din 



(3.4) 



113 



Eq. (|a.4[ ) is the equation for the system of non interacting Pomerons, and the general solution has the 
following form 



N(Y-Y'; 7R ) = (-l) n C nl %G n (Y-0) 



(3.5) 



71=1 



where the coefficients C n can be found from the initial conditions, namely, from the expression for the 
low energy amplitude. For summing the enhance diagrams the initial condition 



GR 



MPSI3 



MPSI4 



N(Y-Y' = 1R ) 



(3.6) 



MPSIIC 
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generates C n 



1 and the solution is 

N(Y-Y'- 1R ) 



1 + lRe A(Y-Y>) 



(3.7) 



31 IC 



The initial condition of Eq. ( |3.6| ) has very simple physics behind it, and has been discussed in Ref. [54]. 

The main idea of the improved MPSI approximation, is to replace 7/j in the generating function 
N ('Jr\Y ) by the low energy amplitude for the dipole-dipole interaction, (see Fig. g). It is easy to see 
that the amplitude in the MPSI approximation has the form (in this equation we denote N (^r\Y) by 
N MFA (7fl|y) where MFA stands for mean field approximation) 



tMPSI 




(3.8) 



MFA 



(2) 



l 7 ( 2 ) =0 7o 



(2) 



where 70 is the low energy amplitude for parton scattering (we will denote this amplitude as 7 and hope 
that it will not cause any inconvenience). 



MPSI5 



3.2 The range of energy for which our approach is reliable 

f ail 
Wf. The most dangerous term is the 

last one. It has an extra Y , which stems from the region of integration yi — 1/2 53 1/A, and it is the first 

term of the renormalization of the Pomeron intercept. One can see (see Refs. [54,55] for details) that the 

renormalized intercept A^j = A — A7 2 *. This cannot be calculated in the MPSI approximation. Therefore, 

the first estimate for the range of energy where we can trust the MPSI approximation, comes from the 

demand that the renormalization of the Pomeron intercept should be small. This leads to 



A 7 y < 1; 



or 



Y < 



1 

A7 



(3.9) 



However, a more restricted region can be obtained from Eq. (p.ib). The term —ju 2 d 2 Z/du 2 describes 



the four Pomeron interaction, and as we have discussed, should be considered as small in our approach. 
Solving this equation without this term, and taking it into account as a perturbation, it is easy to show 
that the contribution of the four Pomeron interaction turns out to be small for 



7F < 1; or 



7 



(3.10) 



*To understand, this it is sufficient to compare this term with the exchange of one Pomeron, which has the form 
51 52 7 exp ( — AV) in our notation. 
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Figure 7: The exact Green's function of the 
Pomeron versus Y — ln(s/s ) for s = 1 GeV 2 , and 
T (Y) for A = 0.339 and 7 = 0.0242. The values of 
parameters have been taken from our fit [7] . 



compar 



Figure 8: The comparison of the exact solution of 
J CFEQ 

Eq. (p?l~b| ) for the Pomeron Green's function G (Y) 
(see Ref. [50,51]) with in the improved MPSI 

approximation (see Eq. ( 3713 )). The figure is taken 
from Ref. [51]. 



The second restriction of our approach, comes from the fact that we consider a'p = 0. Inclusion of a small 
a'jp will destroy our approximation for energies which we can find from the condition 

1 



a'jpY > 1/mf or Y > 
Thus, we can trust the MPSI approximation in the region 



(3.11) 



ER3 



Y < min < — 



1 



1 



7 a'jpmf 



3.3 The Pomeron Green's function and the elastic scattering amplitude 



316 



Using Eq. (ftTqjwe can obtain the Green's function of the Pomeron in a closed form, namely [7] 



G{Y) 



exp 



T (Y) J T(Y) 



r 



T(Y) 



with 



T(Y) 



7 e 



.ay 



(3.12) 



(3.13) 



(3.14) 



and r (0, 1/T) is the incomplete gamma function (see formulae 8.35 in Ref. [58]). Using this function we 
can write the opacities in the two channel formalism (see Eq. ( ]2t|| j in the form 



= 9i(P)g k (b) G (Y) 



ER4 



ESI 



ES11 



(3.15) [ES2 



It should be stressed that the Green's function of Eq. ( |5.13| ) has been found for high energies, where 
the MPSI approach is correct. However, the theory given by the functional of Eq. ( p71 1[) has an analytical 
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solution (see Refs. [50, 51] ) at arbitrary values of the energy. Comparison with the improved MPSI 

f mpar 
shows that this approximation describes the lowest energy within accuracy of 5 

to 10 %, for 7 = 1/16. In our fit [7] 7 = 0.0242 which leads to even better accuracy. Therefore, we can 

safely use the improved MPSI approximation starting from s = 400 GeV 2 which was used in our fit. 



The elastic amplitude is equal 

a e i{b) 



1 [a 



A ltl + 2a 2 (3 2 A 1:2 + /3M 2)2 ) 



(3.16) 



where Aj^ is given by Eq. (JjZ.q ) with fi^ from Eq. ( |3.l5|) . 



■mgrf 



The Pomeron Green's function tends to unity at large values of the argument T. In Fig. |V| we plot this 
function as well as T (Y) for A = 0.339 and 7 = 0.0242 these values were found in our fit (see Ref. [7] and 
below). One can see that G(Y) is quite different from the one Pomeron exchange, even in the region when 
T is smaller than 1. 



hst 



In the MPSI approximation, the exact triple Pomeron vertex (see Fig. y[j 1 is equal to the 'bare' vertex. 
One more comment should be made: Eq. (|3.l3|) leads to 



G(Y) 



00 

£ 

n=l 



(-l) n+1 n!T n (Y) 



(3.17) 



which is a typical Borel summable asymptotical series. Such an amplitude cannot be obtained as a solution 
of the equation after a finite number of integrations. 

3.4 The exact vertex and diffractive production 



ES3 



ES4 






Figure 9: The first diagram for single diffraction. Figure 10: The diagrams that contribute to the 



The wavey lines and the blob denote the exaEdi2 



exact vertex. 



S( jl I Green's function and exact vertex for Pomerons [7] 



The process of single diffraction has not been taken into account in the Green's function of the Pomeron. 
In our paper of Ref. [7] we calculated the first diagram for the single diffraction dissociation, namely one 
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in Fig. ^ In this diagram both the Pomeron Green's function, and the vertex are exact. We found that 
this contribution leads to a small cross section for single diffraction in the region of the large mass, and we 
view this calculation as just the first attempt to solve the problem of the summation of all diagrams. We 
will solve this problem in the next section. The cross section for single diffraction has been calculated in 
Ref. [7] using the MP SI approximation. Here, we calculate the exact vertex shown in Fig. |ftjf We can use 
Eq.3.34 of Ref. [7], where we only need to extract the contribution from the external Green's function of 
the Pomeron, and introduce three rapitities as shown in Fig. |llf The final equation has the form 

r MP5/ (yl,y2,y3) = (3.18) 
x / i \n+m an am 

V ^— -y n + m — - V MFA (w P W P -vl-V )\ 1-1 

n=l:m=l 



d n V MFA ( w t ) y m _ y2 ) i d m T MFA (w\y m -y z ) 

Qn w t 

A 1 



Qn w t w 1 d n W t W 1 



where 



and 



, r , . - r( r {r 1 (2T(y 1 -y 2 )) - T x (2T ( yi - y 3 ))} (3.19) 

4 T(yx - y 2 ) - T(y 1 - y 3 ) 

with Ti (T) = (1/T 3 ) x {T(l + T) - exp (-1/T) (1 + 2T) T (0, l/T)} (3.20) 



T MFA {w,w; yi -y m ) = 2Aww- (3.21) [smT 

(1 + (w + w)(e A ^-y^ - 1)) 

rMFA <» ; y - ») = + p- 22) ^ 



where w and w are the variables that we needed to introduce, as has been explained in Ref. [7]. It is easy 
to see that 

T MPSI (yl,y2,y3) ^Z^iBZ^ ^ e -A(2 yi -y 2 - ys ) (3>23) 

and therefore, to a good accuracy we can consider that the exact vertex is equal to the bare one, namely, 
r MPSI (yl,y2,y3) = A. 

The double diffraction contained in the Pomeron Green's function has been calculated directly using 
the unitarity constraint (see Ref. [7]). 



SD5 



4. Summing the full set of diagrams 
4.1 Elastic scattering 

In this section we sum the full set of the diagrams using the improved MPSI approximation. In our 
approach that has been discussed in the previous section, the elastic amplitude is written as the sum of 
the eikonal diagrams for each state i in the two channel model, with the exact Pomeron Green's function 
(see Fig. [1 i\- A) . Here we sum all diagrams (see Fig. |11|-B using the MPSI approximation. The diagrams of 
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v 



f ulset 



Figure 11: The full set of the diagrams. Fig. ffr The sum of enhanced diagrams in the two channel approach, 
Fig. plpB shows the full set of the diagrams which in Fig. pTpC is pictured in a way that is most suitable to illustrate 
the MPSI approach. The bold wavey line stands for the exact Pomeron Green's function that includes all enhanced 
diagrams. 

Fig. |ll|-B (in the MPSI ap prox imation) , are shown in Fig. |ll|-C. The first step in determining this sum is 
to find the solution to Eq. (p.4|) with the initial condition 



set 



Jv in ( 7; y-y' = o) = i -e 



-9i7 



(4.1) 



ESEl 



The meaning of this conditions is very simple and stems from Eq. (|2T9| ) : Eq^ £T f gives the eikonal formula 
for the scattering amplitude for Y — Y' = 0. The general solution of Eq. (|5t!)is simple, namely, 



N 



MFA 



(Y-Y'; 7R ) = N in ( 1 = N(Y-Y'; lR )) 



(4.2) 



3J4 



where N (Y — Y'\ 7^) is given by Eq. (|3.b[ ). 
For the initial condition of Eq. ( j^lj ) 

™ FA {Y-Y'- 1R ) = 1 -exp{-ffi 



N 1 



1R 



e A(Y-Y>) 



1 + lR e A ( Y ~ Y ') j 

Using the generating function for Laguerre polynomials (see Ref. [58] formula 8.973(1)), namely 



(4.3) 



(1-z) 



-a-l 



exp 



X z 
z - 1 



(4.4) 



ra=0 



FSE1 



FSE2 



FSE3 



LPGF 



we obtain for Eq. (gTa) 
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N MFA {y _ y,. lR) = _ £ L -l Cgi) ( lRe ^Y-Y>)Y (4.5) 



FSE4 



n=l 



Using Eq. ( |I ? b[ ) and Eq. (^^pwe have for the scattering amplitude 



iV™(Y;) = ^n\L^(~ gi ) L^(g k )(- 7 e AY ) n (4.6) [fses] 

71=0 

Introducing n\ = j Q °° £ n d£ we can re-write Eq. ( jljfij ) in the form 

JV MW (F;) = / ^e-^^ L" 1 ^) L" 1 (£ fe ) (-^oe Ay )" (4.7) [FSE6] 

n=0 

Using formula 8.976(1) of Ref. [58], namely 



n=0 

we derive the final result 



SUML 



N - k (y) = J 7 i + mr) exp i" eT(y) rrw) / Jl ( 2 ttwt J (49) ™ 

where 

T(Y) = 7 e Ay where 7 = 7o . (4.10) 

Eq. (fhsjj) reduces to a simple and elegant formula in the case that we can consider giT(Y) ~ 1, 
gigiT{Y) > 1 but T(Y) <C 1. Indeed, in this case the integral over £ can be taken and 

A ik (Y;b) = 1 - exp(- i;t fcr r (y) -\ (4.11) H^o] 

' v 7 I 21 + T(y) & + &] J 1 ; 1 1 

In Eq. ( §§ T we neglected 6 dependence. In our approach, the b dependence enters through the vertices 
gi(b). For Eq. ( ^.11| ) it is easy to write the expression that takes into account the correct impact parameter 
behaviour, namely, it has the form 



/ d 2 */ V W r ) '—. '—^ ) (4.12) [FSE91 

j 1 + rm o, 1 v 1 -+- m. 1 b> 1 



/c (F; b) = 1 — exp < 

1 + T(Y) [& (yj + g k [b- V 
To obtain the elastic amplitude we need to substitute the amplitudes in Eq. ( |lllG| ) 
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simsetim 



Figure 13: The set of diagrams that is selected and summed using the fact that giG (T (Y)) > 1 while 
A 2 G(T(Y)) « 1. 




For the nucleus- nucleus scattering Eq. ( pLTi ) is the gen- 
eralization of the Glauber formula, and we intend publishing 
our estimates for the cross section of nucleus-nucleus scatter- 
ing elsewhere. 

FSF.ft , 

Eq. (gTTTj) is very simple and we intend to use this ap- 



proximation for the description of the experimental data. 

icmiparl 

We plot in Fig. f2 the comparison of ImA e i (Y) given by 
the correct Eq. (£13 J" , with the approximate Eq. ( |4 for 



7 = 0.0242 and A = 0.339 at the values of g\ = c/2 = 3 and 5. 
The values of parameters 7 and A are taken from our fit [7] 
with th e amplitude that takes into account the emhanced 



C( ?fiS^Fai Qs only (see section 3). Note that at large g, the ac- 



Figure 12: Comparison of the exact imagi nary curacy of the approximate solution is about 5%, while at low 

part of the elastic amplitude, given by Eq. (pH^ ' <•-.-. n,, , - n n-, , ^. l^ffP arl 

with the approximation for gT(Y) « 1 while 
AT (Y) < 1. The values of 7 = 0.0242 and 

A = 0.339 were taken from the fit [7] using the Pomeron diagrams, without assuming the MPSI approxima- 

sum of enhanced diagrams (see section 3) tion (see Ref. [8]). 



values of g the errors could be as large as 12%. (see Fig. 

Eq. (gTTJj can be derived by direct summation of the 



Using the fact that Eq. sums the net set of the 



diagrams not only in the MPSI approximation, we can also sum more complicated diagrams in which the 

IfjjJ set JESj— . teigsgtim 
'bare' Pomerons in Fig. |11|-C are replaced by the exact Pomeron of Eq. Q3.13D (see Fig. |13|, where we show 

examples of the diagrams that we sum, as well as examples of the diagrams that we still need to calculate. 

In the selection of these diagrams we used parameters: gi^G (T (Y)) > 1 and A 2 G(T(Y)) <C 1. In the 

fit of Ref. [7] g > 5.8 and A = 0.339 and, therefore, we can obtain the scattering amplitude using this 

re-summation procedure within the accuracy A 2 /^ < 0.1/5.8 = 0.02. We need only to replace T (Y) in 

Eq. (gf§ by G(T(Y)). 



4.2 Diffractive production 

To calculate the single diffraction process it is necessary to introduce three different variables : 7i?,7i? and 
"jin to describe Pomerons in the amplitude, in the complex conjugate amplitude, and the amplitude of 
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the cut Pomeron, respectively (see Ref. [54] for details). The difference to the calculation of theelastic 
amplitude lies mostly in the fact that we need to take into account all three amplitudes in Eq. ( fT^t ). For 
the case of single diffraction, we have only one cut Pomeron at Y = Y m that decays in one 7^ Pomeron, 
and one 7i? Pomeron. Therefore, the cut Pomeron in Fig. Jl^fhas the following form [54] 

«Pomeron(^^;F-F M ^y m = ln(MV S o)) = 2A TO 7k (i - j£^^f ^ ^ 



JR 



iR 



Y' 




1R 



Figure 14: The MPSI approximation 
for the cross section for single diffrac- 
tive production of mass (M 2 ,Y — Ym = 
ln(M 2 1 so)) . The dashed lines show the 
cut Pomerons. All other notations, are 
as in Fig. IT. 
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Choosing Y' = Y M we have the N MFA (jr,j r ;Y - Y M = Y m = ln(M 2 /s )) for the full set of the 
diagrams in the form 

N MFA ( lR ,j R ;Y-Y M = ln(M 2 /s )=Y m ) = (4.14) 

N cut Pomeron fa. ™ Y ~ Y M = m(M /s ) = Y m ) exp (-— — - ^- j exp — — j 
Using Eq. (0713 ^ we can reduce Eq. ( |4,l4[ ) to the form 

N mfa {ir,7r\Y-Y m = HM 2 /so) = Y m ) = (4.15) 

2Ajp7r7 R tdte y ^ R ^ R > p ^ m » exp -— -^r- exp - _ — ^y- 

Using Eq. (jojfand Eq. ([hbjj we can expand N MFA with respect to powers of 7^ exp (AY m ) and y R exp (AY m ), 
namely 



N MFA ( lR ^ R ;Y-Y M = ln(M 2 /s ) = Y m ) = (4.16) 

j.k+k 



2A* EEEffl ^* (- 4y -)" ^V. (-5.) ("ft) 



jV=n+n fc=0 fc=0 

where n, k, n and are integer numbers. 



- 19 - 



Using formulae 8.972(1) and 9.211(2) of Ref. [58] we obtain 

J ESSD5 J FSSn .4 

Substituting Eq. ( H7T7 ) in Eq. p.lq ) we have 



t k+k , „ a , ak,\JV 



iV MFA (7K,7^ m ) = 2Ap,e-^ £ £ £ ^ (-1)" 7 " ^ (-e AY ™) 

N=n+n k=0 k=0 

1 



(n-k- 1)! 
1 



-9i 



(n-k- 1)! 

The generating function that describes the low cascade is equal to 

N MFA (rf R , jr; Y — Y m ) = N MFA ( 1R] Y - Y m \Eq. (ff) x N MFA ( 7 r; Y - Y m \Eq. (gf) 
and can be re-written in the following form 

n N 

- E Ei- 1 )" h A{Y - Ym) ) 

N=n+n 



MFA 



N 



1 

x — e 
n! 



■;?■! 



(-<7k) 



^"e-«r +f Ji ( 2a/ M 



iI6 



Eq. ( |3.8| ) should be replaced by a more general equation, namely. 

("l) n 



AT MPSI 



(Y m , Y Y m 



n=l 



n! 



E (7o Tsd)' 



n=l 




MFA (y . (1) -(1) 
I m,lR ,1r 



MFA 



Y -Y r 



■ J 2 ) ^,( 2 ) 



Using Eq. ( fi.lYf) we can sum over A; and k reducing Eq. (|4.1q ) to the following equation 



N MFA U£^;Y m 



2A^^e-^T§ D ~ gi e-^ £ 



-1)" ^ (_i)f 



E 



1 ("-1)1 ^ (n-1) 



(4.17) 



FSSD5 



(4.18) 



(4.19) 



(4.20) 



FSSD7 



(4.21) 



(4.22) 
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Using Eq. fl4.22p we can reduce Eq. to the form 

N™ D PSI (Y m ;Y -Y m ) = 2A~g i7 2 e- AY -T 2 SD (Y,Y m ) J tdte'^^ J d& d& -^L= (4.23) 

x E (^rw (* + ^) n ~ x £ _1 Ji Ji (2v^6) x {6 6; & -> ^} x (- 7o T sl) ) n +" 

n=l 



Summing over n and n we obtain 



ATjif S/ (^m;^-^m) = (4.24) 
2A e - Aym Tj D (Y,Y m ) g-i e^ 7 *'- 2 ** j tdt J d£ 1 d&dtid&e- t & + ^ TsD - t e-t 1 -t 2 -t 1 -b 



and after interation over t we have 

A^'O^Y- Y m ) = (4.25) 

2A e- AK - T| D (Y, Y m ) ~g^ e~ 2 '^ +2 ^ [ d ^ d ^ ^ e -fc (4.26) 

•/ (1 + (6 + 6)7 t 5d) 

x -^-Le-^^fe^e^^ 2 ^ J: (2v^6) Jx (^y^T) J x (2^2") Ji (2^2) 

Introducing new variables y = £1^2 and y = £1^2 and integrating over them, we obtain the final formula 
for the single diffraction amplitude Aik (Y, Y TO ) in the form 

A™ (Y, Y m ) = 2A e~ AY - T 2 SD (Y, Y m ) ~ 9l v 7/^ 2 * * " [ ff»«fey^6 -6-6 (4 . 27) 
x * (2V&&) J, (2^) {l - e X p (n^-) } {l " exp ( TT |y } 

where 

T 5D (Y;Y m ) = 7 (e AYm - 1) e A ( y - y ™) (4.28) [tsd] 



In Eq. (gT27j) we took into account that in the two channel model, the low cascade could be initiated by 
different states (k and I ). As in the case of the elastic amplitude, gi in Eq. (pL2Y| ) should be replaced by 
gi (b'^J and by g^^ (b — b'j and should be integrated over d 2 b' . 

The total cross section of the diffractive production can be wriiten as a sum of two terms: the Good 
-Walker term which is equal to 

a GW = j d % ( a p{- a 2 Af A + (a 2 - /3 2 ) Af >2 + /3 2 A% 2 } ) 2 (4.29) | FSSDGW 
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where A^k are given by Eq. fl4.1l| ), and the term which describes the diffractive production in the region 
of large mass, namely 

^Large mass = 2 J dY m f d 2 b (4.30) 

{ a 6 Afg A e~ n ^ Y ^ + a 2 /3 4 4? i2 e ~ n ^ Y ^ + 2 a 4 /3 2 Af° 2 e -U n ^)+^(y,b)) 
+ /3 2 a 4 Al° ie - n ^ ¥ V + 2/3V A SD 2e -l(nMY;b) + n 22( Y;b)) + ^ ^ e -n 22 {Y-,b) j 

To find the cross section for double diffraction we use the s-channel unitarity constraints as was 
suggested in Ref. [7] 

2Af !k (Y;b) = \Af !k (Y;b)\ 2 +2A?£ k (Y;b) + A^ + A%{Y;b) (4.31) [755511 

It has been shown in Refs. [47,54,59] that the inelastic amplitude Af^ k (Y;b) = Af k (2T(Y);b). Therefore, 
from Eq. (tHjjfr the amplitude for double diffraction production is equal to 



A° k D (Y;b) = (4.32) 
2Af >k (T(Y);b) ~ \A^(T(Y);b)\ 2 -2 / dF m Ag£ fc (T SD (Y, Y m ; b) ; b) - A\\{2T (Y) ;b) 

Finally, the cross section of the double diffractive production is the sum of the Good- Walker contribution, 
which has the form 

<r$g = J d 2 b a 2 (3 2 {Af tl - 2Af >2 + A^ (4.33) 

with Af k given by Eq. ( |4.1l| ), and the term which is determined by the Pomeron interaction, and which 
contributes to the production of large masses, namely, 



a 



Large mass 



DD 



5. Corrections to our approach 

5.1 a'jp ^ and the Pomeron as a fixed branch point 

In constructing our model we have made two assumptions that considerably simplify our approach. Firstly, 
we considered the case of a'jp = for summation of the enhanced diagrams, in spite of the fact that 
a'jp ~ 0.02 GeV ~ 2 , is the value obtained from our fit of the experimental data. Secondly, we replaced 
the Pomeron Regge pole by a fixed square root singularity, that appears in the N=4 SYM. Both these 
assumptions are not principle in nature, but have been made to facilitate the simplicity and transparency 
of our approach. In this section we will expand our formalism so as to assess the consequence of our 
assumptions. 



FSDD3 



J d 2 b {a 4 Aff e ~ n ^ Y ^ + 2a 2 /3 2 Ajf, e - Q ^ Y ' b \ + /3 4 A°g e~^ Y ^ } (4.34) [fsdm] 
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/ a / a We know how to include the impact parameter^ 

dependence in the generating function formalism, Eq. ( j^jj ) 
should be replaced by a new equation [48], which has 
the following form 

(0,0) £ , "2, (o, o) . (y-lo;6;u) 




eqpea- 



<9F 



(5.1) 

Figure 15: The graphic form of the equation that ^ v ^" (y _ y fc } _ Au (i _ u) Z (y _ y o; . 
sums 'fan' diagrams. The Green's function of the ex- 
act Pomeron is depicted as bold wavey line, while the JCQE1 

'bare' Pomeron is shown by thin wavey line. E Q- (CT " is written for the generating function that 

sums the 'fan' diagrams. For the scattering amplitude, 
this sum is the solution of the simple equation shown in Fig. |lo|, and it is of the form 

G (y; b) = Go (y, b) - G 3]P £ dy' j d%' G (y - y',b - fj G 2 (y 1 - 0, b') . (5.2) 

Go (y, b) is the Green's function of the 'bare' Pomeron 

G„M) = -A-*'^ ^ * i2) (b), (5-3) 
Aitotjpy 

which is the Fourier transform of G (y, q) = exp (Ay — a'jp g 2 ) , where q 2 = —t is the momentum transfer 
squared. Go satisfies the equation 

dG °Q V y h) ~ oc'jp V 2 Go (y, b) = A G (y, b) . (5.4) 

Using Eq. (irjf) , Eq. can be rewritten in the form 

dG{ Q y h) ~ a! P VlG{y,b) = AG(y,b) - AG 2 (y,b) , (5.5) 

where we substitute G^p = A as in Eq. 

gfjRl 

This equation is the same as Eq. ( |b.l[ ) for the generating function. Therefore, in the framework of the 
MPSI approximation Eq. (p together with the obvious property of 

Go (y, b) = j d 2 b' Go (y -y',b- fj G (V, V) , (5.6) 

lead to the Green's function of the exact Pomeron that includes the impact parameter dependence. 

As we have discussed, the solution to Eq. ( j^ffi can be found in the form Z (Go (y) f(u)) in which f(u) 
that follows from the equation, and the form of Z stems from the initial condition Z (y = 0, u) = u. The 
solution is 

Z(y,u) = — —- . r . (5.7) 

u + (l-u)Go(y) 



CDR2 



CDR3 



CDR4 



CDR6 



CDR7 
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|f nRi |f nRj 

We suggest that the solution of Eq. (p.lp has the same form as Eq. Qb-Vp , namely, 



Z(y,u;b) 



u 



u + (l-u)G (y;b)' 



(5. 



One can see that this form of Z(y,u;b) satisfies the initial condition Z (y = 0,u;b) = u5^ [b] and 
Z (y,u = 1; 6) = 1. Inserting Eq. (Ib^) one sees that it does not satisfy Eq. ( fxlj) . We consider 
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dZ (Y-Y -b-u) 
dY 



(5.9) 



a'jpV 2 b Z (Y-Y ;b;u) + Au(l-u) Z (Y-Y ;b;u) 



Gl (y,b, 



d 2 Z(G (y,b);u) 



x 



/2 „,2 • 



^From Eq. ( p.9| ) we see that if b 2 < (a'p y) , the cor- 
rections are small. Recall that the Pomeron contribution 
is dominant at b 2 < Aa'p y. It is also small at large values 
of y since Go S> 1. 



23456789 10 

b in 1/GeV 



Having Eq. ( |h.«|) as the solution for the generating 
function, we can find the exact Green function (G(y,b)) 
aitgi the Pomeron in the MPSI approximation which is 

Figure 16: The scattering asmplitude A n and ec l ual to Ec h ( E^D witn 
A e i for a'p = (solid lines) and ap = 

0.01 1/GeV 2 , (dashed lines). All other parameters T <T) — ► T ( Y 'i b ) = ( 5 - 10 ) 

are taken from Rcf. [71. 1 Ay 

= 7G0M,) = 



4 a' j/ 



Instead of Eq. ( 



Ina'jpy 
for fijfc we have 



In Fig. ^we plot the b space amplitudes for the highesr LHC energy with a' F = and a^, = 0.01 GeV~ 2 . 
The difference is negligibly small and, therefore, we neglect the restriction on the kinematic region that 
followed from taking a'p = (see Eq. ( flRiP and Eq. (pTYj )). Note, our formalism ( due to the MPSI 
approximation) is only valid for W < 100 TeV [7], this is far beyond the LHC energy range. 



Eq. (p. 4 ) shows that a method to include the energy behaviour of the 'bare' Pomeron, is different 
from the one for a Regge pole. As we have mentioned in N=4 SYM the leading singularity in the angular 
momentum plane is not a pole but a fixed cut, with an energy behaviour 
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(5.12) C0R12 
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Figure 17: Contributions of secondary Reggeons denoted by zigzag lines. Wavey lines denote the Pomeron. Fig. [Hp 
A shows the contribution to the scattering amplitude due to exchange of a Pomeron and Reggeon. Fig. [lYkb and 
Fig. |i y[-C illustrates the fact that the exchange of secondary Reggeons can be reduced to the inclusion of the vertices 
for Pomeron-Pomeron interactions, or can be included in the Good- Walker mechanism. 



which satisfies the equation 

dG (y,b) 
dy 



AG (y,b) - j /yG (y,b) -> AG (y,b) 



(5.13) 



^From Eq. (p.lj ) one ^an see that we can neglect tie contribution of the factor, y J ° ri by o differentiating 



over y. Therefore, Eq. ( |o.2p can be reduced to Eq. (|b.a| ) which has solution of Eq. flb.Sp with Go from 
Eq. (CT - 
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5.2 Low energy description and the 'threshold effect'. 



5.2.1 Low energy behaviour of the scattering amplitude 

As we shall see in the next section, most of our data base consists of lower energy points from ISR and 

SppS/SppS (W w 20 -70 GeV), where the contribution of the secondary Regge exchanges are important. 

A secondary Reggeon has an energy behaviour exp (Ajj(Y — 0)). The sum iP + M describes the energy 

behaviour of the elastic scattering amplitude without screening corrections. This sum replaces the single 

Pomeron exchange in the definition of Q^. Inserting this sum everywhere in the more complicated diagram 

(see Fig. LV[), one can see that the integrations over rapidities reduces the contributions of the secondary 

— ■ • n "pp? 1 

Reggeons. By introducing new vertices for the Pomeron-Pomeron interactions ( see Fig. ™ the inte gra- 

tion over y\ — yi can be replaced by a new P — > 3JP vertex), or it can be absorbed into G-W mechanism 

iXoaen 

(see Fig. |1V[-C Since in our approach the vertices, other than the triple Pomeron vertex are considered to 
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be small, we arrive at the conclusion that for lower energies we only need to replace the single Pomeron 
exchange by IP + M, in the definition of Ojfc. 

[settle 

In Fig. |18| we compare our prediction for 
lower energies with the experimental data and 
obtain a satisfactory description to within 10%. 
The conclusion is very simple: we do not need 
an additional source to describe the lower energy 
behaviour of the amplitude. 

5.2.2 The 'threshold' effect 

As a consequence of the above conclusion there is 
no requirement to introduce a Pomeron thresh- 
old, i.e. we do not need to assume that the 
Pomeron contributes to the scattering amplitude 
for Y > yo with yo sj 1.5 -f- 2.5. 
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Figure 18: The total cross section (atot = 
1/2 (crtotipp) + ctotipp))) a t l° w energies in the frame- 
work of our approach with the parameters determined by 
high energy data. The curve illustrates our parametrization. 



influence of the threshold cutoff on the value of 
the scattering amplitude, and especially on the 
value of the survival probability. We calculate 
the simplest enhanced diagram of Fig. 



ducing this cutoff, assuming that Y 



well as U2 > yo- We obtain 



A(Fig.BJ 



ffV -2A(Y- 

A 2 



■yo) 



-AY\ 



mtro- 

yi > y as 

(5.14) 



First, we would like to draw the readers attention to the fact, that without a cutoff the typical Y — y\ ~ 1/A 
and y2 ~ 1/A. Therefore, we need to introduce a cutoff only if yo > 1/A. Second, we need to change 
7—7-7 exp (A yo) and multiply the Pomeron exchange by G (Y — yo)- Doing so in our parametrization we 
obtain the following values for (\S^ nh \} for yo = 0,1.5,2.3. For the Tevatron energy (W = 1.8 TeV) we 
obtain 0.285, 0.7 ,0.99. For LHC (W = UTeV) the corresponding values for (\S^ nh \) are 0.06, 0.12, 0.19. 
(|$? n J) is the survival probability initiated by the Pomeron interactions. The conclusion from this exercise 
is very simple: in our approach we do not need to introduce a threshold for Pomeron exchange. 
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6. The fit to the data and its phenomenology 



6.1 The main formulae of the fit 

The main formulae, that we use, have been given in sections 4.1 (Eq. ([4.12 )" and 4.2 (Eq. ( |4.2Y ),Eq. (gTT^jTEq. (| 
and Eq. ( J0$4| j). However, we change these formulae so as to take into account the enhanced diagrams. As 
has been discussed in section 4.1 this procedure reduces to the following substitution in our basic formulae: 

T(Y) (Eq. ffik) — >• G(T(Y))(Eq. 



T SD (Y;Y m ) (Eg. 



— G(T (Y n 

7 



); 



G[T (Y — Y m ) ; 



(6.1) 



We also prefer to use gi and G^jp, instead of gi and 7 in the main formulae of Eq. ( |4.12 f, considering 
only GIT (Y) I as a function of 7. Finally, this formula has the form 



Ai, k (Y; b) 



1 



^(Y;b) 



d 2 b' 



1 — exp 
9i(b') g k [b-b' 



2 «F C^; b) 



1/ 7 G(T(T)) 



1 + {G$p/i) G[T(Y) 



9i [V] +9k [b-b 



(6.2) 
(6.3) 



Note that gi in Eq. flb.ap have dimension of inverse momentum (see Eq. ( p.!9| )) as well as G^jp, while 
7 is dimensionless. Actually 7 2 = J* d 2 kG 2 jp, but because we do not know the dependence of G3JP with 
respect to transverse momenta of Pomerons, we consider 7 and G^jp as independent parameters of the fit. 

In Eq. ( jdlz^D we have to replace all gi by and multiply it by factor G3JP/7 2 in addition to the 



substitution of Eq. 

For completeness of presentation we list below the formulae for physical observables (see Re. [7] for 
details). The amplitudes for the observable processes have the form 

a e i(s, b) = i{a 4 A hl + 2a 2 (3 2 A^ 2 + /3i4 2 , 2 }, (6.4) 

a sd (s, b) = iaP{-a 2 A lA + (a 2 - f)A l>2 + p 2 A 2y2 }, (6.5) 

a dd = ia 2 p 2 {A lA - 2Ai 2 + A 2j2 }. (6.6) 



FIO 



FIMF 



EL 
SD 
DD 



The corresponding cross sections are given by 

o to t{s) = 2 

CTsd(s) = 



d 2 ba e i (s,b) , 
d 2 b\a e i (s,b) | 2 , 
d 2 b\a sd (s,b) | 2 , 
d 2 b\a dd (s,b) | 2 . 



(6.7) 
(6.8) 
(6.9) 
(6.10) 
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XSEL 



XSSD 



XSDD 
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6.2 The strategy of our fitting procedure 



Our Pomeron model, as well as the KMR [9] and Ostapchenko [10] models have the same following ingre- 
dients: 

1) A bare non-screened Pomeron exchange amplitude. 

2) s-channel unitarity is enforced by eikonal rescatterings of the colliding projectiles. 

3) These rescatterings proceed through elastic and diffractive states according to the G-W mechanism. 

4) t-channel unitarity is maintained through the Pomeron interactions. 

In our model the Pomeron is specified by nine parameters, the Tevatron data on its own is not sufficient 
to determine the parameters. Consequently, we have also to include ISR - SPpS/SppS lower energy (W 
~ 20-70 GeV) data. This data has small errors which facilitate a reasonably reliable fit. To reduce the 
number of Reggeon parameters, we define a tot = \{&tot{pp) + &tot(pP)- The inclusion of the Regge sector 
of our fit requires five additional parameters, i.e. we have fourteen parameters in all. 

Our data base has 58 experimental data points, which include the p-p and p-p total cross sections, 
integrated elastic cross sections, integrated single and double diffraction cross sections, B e i, to which we 
have added a consistency check of the CDF data (-t < 0.5GeV 2 ), d ^ff /s ( t = 0.05 GeV 2 ) and B sd . 

The data points were fitted to determine the 14 free parameters of our model. We fit simultaneously 
the entire data base. The only minor tuning which we employ is discussed in section 6.2. 

As was mentioned in the previous section, most of the experimental data is available at low energies 
[yfs ~ 20 — 70 GeV) where the secondary Reggeon contributions are essential, this data has small errors. 
We deal with the secondary reggeon in the same way as in Ref. [7] adding to fijp in Eq. ^>.% . 

Using the parameters of the fit (see Table 1), we find that the contributions of the large mass diffraction 
to the single diffractive cross section as well as to the double diffractive cross section, are rather small. 
Therefore, we neglect these contributions and use Eq. ( [4.29 pand Eq. (f03j in the fit. After determining 



the parameters of the fit we use them to describe the large mass diffraction. 
6.3 The results of the fit 

As stated above, our fit is based on 58 experimental data points. The model gives a good reproduction of 
the data, with a x 2 /d.o.f. = 1.56. However, a large contribution to the value of x 2 /d-o.f. stems from the 
uncertainty of the value of two single diffraction cross sections, and of the CDF total cross section [57] at 
the Tevatron (W = 1800 GeV ). Neglecting the contribution of these three points to the total x 2 we obtain 
X 2 / d.o.f. = 0.86. The quality of the description of the experimental data is demonstrated in Fig. |itf The 
values of fitted parameters are listed in Table 1. An important advantage of our fit is that it provides a 
good reproduction of add- 



In Fig. 2C the amplitudes A^. are plotted. One can see that in spite of the smallness of a'jp we reproduce 
the growth of the radius of interaction with energy. The values of physical observables for higher energies 
are shown in Table 2. 
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Figure 19: Comparison w_ith the experimental data the energy behaviour of the total (Fig. |l!d|-a), elastic^ ^Fig. |l!frb), 
single diffraction (Fig. |ft)|-c) and double diffraction (Fig. ^|-d) cross sections and elastic slope( Fig. |lu[-e) . The 
solid lines show this fit while the dashed lines correspond to the fit of Ref . [7] . Fig. |lt|-f shows the behaviour of the 
amplitude Ai^ as function of the impact parameter b for the Tevatron energy. 



6.4 Comments on the parameter values of the fit 

An attractive feature of our fit is that the exceedingly small value of a'jp found in the original fit [7] is 
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a'jp 


9i 


92 


mi 


m 2 


0.2 


0.388 


0.020 GeV~ 2 


2.53 GeV- 1 


88.4 Ge^" 1 


2.648 GeF 


1.37 Ge^ 


A K 


7 


a 'm 


9? 




^0,1 


G3JP 


-0.466 


0.0033 


0.4 GeV- 2 


14.5 GeF- 1 


1343 GeV' 1 


4.0 GeV- 2 


0.0173GCF- 1 



[tl] 



Table 1: Fitted parameters for our model. The quality of the fit is x 2 jd.o.f. = 0.86 (see the detailed explanation 
in the text) 



reproduced. However, the values obtained for Ajp and 7 are smaller than our previous values [7]). Our 
results suggest that the complete summation of the Pomeron interaction sector, presented in this paper, 
induces a weaker screeing that was found by the partial summation presented in Ref. [7]. A consequence 
of this feature is that S 2 nh calculated in this paper is expected to be larger than the corresponding values 
obtained in [7] . 

The small value obtained for 7 = 0.0033 is encouraging, since 7 oc a 2 in QCD, supports our key 
supposition that rather short distances, contribute to the soft interaction at high energy, in agreement 
with all approaches to the Pomeron structure considered above. Note, that since 7 2 = J* d 2 kG 2 jp we can 
evaluate the value of the typical transverse momentum of the Pomeron in the triple Pomeron vertex, which 
turns out to be « lGeV. 

g<i is rather large, as in our previous approaches (see Ref. [7]). The consequence of ^ 3> 1 is seen 
in Fig. [l^j-f, where the amplitude are shown at the Tevatron energy. One can see that amplitude 
A22 is equal to one in a wide region of b and, therefore, corresponds to black disc scattering. However, its 
contribution is proportional to /3 4 in elastic amplitude and to (3 2 in single and double diffraction amplitudes, 
consequently, it's relative contribution is small. 

It should be stressed that the value of our phenomenological parameters (see Table 1) are in agreement 
with the theoretical estimates of Eq. ( pTSp . Choosing the typical soft scale [i = lGeV we can see that 
gifj, « 1 and 

G 3iP/ u rj 7 « Af, « Ajp < gm (6.11) pFI3 



As we have discussed, in this paper we sum all diagrams in an approximation in which giG[T (Y) ) > 1 
while A^G^T (Y) ^ -C 1. The values of the fit parameters support the use of the approximation. 
In Fig. jltj the amplitudes are plotted. 

The values of physical observables for higher energies are shown in Table 2. 



6.5 Comparison with other approaches. 

At present there are three groups that are working on modelling the strong interaction at high energies: 
Durham group [9], Ostapchenko [10] and our group. The models have a lot in common: a rather large 
value of Ajp, small values of a'jp, a large contribution of the Good- Walker mechanism, and a significant 
Pomeron-Pomeron interaction. The theories mainly differ in the way the Pomeron-Pomeron interaction 
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Figure 20: The amplitudes versus b at LHC energies. 
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91.3 
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10.2 


6.46 


19.3 


10 


96 


24.6 


10.5 


6.48 


19.9 


14 
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26.1 


10.8 


6.5 


20.5 


100 


128 


35.6 


12.7 


7.79 


29.9 



t2 



Table 2: Predictions for energies that will be accessible at the LHC. 



is taken into account. In the Durham model as well as in the Ostapchenko one, all vertices of Pomeron- 
Pomeron interactions are included, using a slightly different phenomenological anzats. In our model we 
include only the triple Pomeron vertex, as has been discussed in section 1. In the Durham and in our 
models Ajp « 0.240.3 and a'jp « 0, while the Ostapchenko model gives Ajp = 0.14 and a'jp = Q.l&GeV . 
All models give more or less the same results at the Tevatron energy, and their results at the highest LHC 
energy: span the values a tot = 86 4- HAmb,a e i = 20 4 30 mb, a s d = 10 4- 16m6 and add = 54 13m6. 

Note that the Ostapchenko model also includes a hard Pomeron, while our model has only a single 
Pomeron. The KMR model is more complicated as the Pomeron dependence on k\ is described by three 
Pomerons. As one can see from Table 3, the Ostapcheko model leads to the largest value for the total and 
elastic cross sections, and gives small values for the double diffractive cross section. The Durham model 
predicts the lowest value of the total cross section but a rather large double diffraction cross section. 
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8.87 10.9 
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13.3 


16.1 


11.0 


(7 dd (mb) 


5.36 


3.53 7.2 


13.3 


3.93 


6.1 


6.5 


13.4 


12.9 


4.83 



Table 3: Comparison with the other models: GLMM is our model in which we summed only enhanced diagrams [7], 
GLM is the model described in this paper, KMR(07) and KMR(10) arc two models of Durham group ((Ref. [9]) and 
OS(C) is the model developed in Ref. [10]. The predictions of the fit of 2010 is preliminary and are taken from the 
talk of A.Martin at DiffractionTO. 



7. Conclusions 

In this paper we have constructed a model to describe strong interactions at high energy, based on two main 
theoretical criteria: it should include the main result of N=4 SYM, which is the only theory that is able to 
deal with a large coupling constant; and it should provide the natural matching with high enerrgy QCD. In 
accord with these ideas we assumed that Ajp is relatively large and a'jp —> 0. Using these assumptions in 
this paper we sum all enhanced and semi-enhanced diagrams. The enhanced diagrams have been calculated 
in our previous paper [7]. For the first time we obtain analytical formulae for the scattering amplitude 
both for elastic scattering and for diffractive production. 

Using these formulae we made a fit to the available experimental data and predict the main soft 
observables at LHC energies (see Tables 2 and 3). 

We study the accuracy of our approach related to the corrections induced by small (both not equal to 
zero a'jp) and by the fact that the 'bare' Pomeron cannot be a Regge pole, but a cut. 

In this paper we have completed the formulation of a theoretical self consistent approach. In the future 
we intend to apply this approach to the numerous practical problems, such as survival probability for di-jet 
and Higgs production, inclusive cross sections, rapidity and multiplicity correlations. 
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